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Abstract 

We investigate nuclear matter, i.e. the nuclear equation-of-state (EOS) as well as the relativistic 
mean fields in the chiral limit. The investigations are based on a chiral nucleon-nucleon EFT 
interaction where the explicit and implicit pion mass dependence is known up to next-to-leading 
order. The nuclear bulk properties are found to remain fairly stable in the chiral limit. Based on the 
same interaction the in-medium scalar condensate is derived, both in Hartree-Fock approximation 
as well as from the Brueckner G-matrix, making thereby use of the Hellman-Feynman theorem. 
Short distance physics which determines the reduction of the in-medium nucleon mass is found to 
play only a minor role for the reduction of the chiral condensate. 
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I. INTRODUCTION 



The investigation of dense nuclear matter is one of the most exciting topics of present days 
nuclear and hadron physics. Since hadrons are excitations of the QCD vacuum their in- 
medium excitation spectrum is closely connected to the modification of the QCD vacuum 
at finite density. The chiral condensate (qq) is expected to be reduced at finite density 
and/or temperature which leads to a partial restoration of chiral symmetry of QCD and 
should be reflected, e.g., in shifts of the corresponding hadron masses. Chiral symmetry is 
an exact symmetry of QCD in the limit of massless quarks. Since the up and down current 
quark masses are small, i.e. of the order of 5-10 MeV, this symmetry is still approximately 
fulfilled. In nature it is, however, spontaneously broken by the non-vanishing - and large - 
expectation value of the scalar quark condensate (qq) of the QCD vacuum. The spontaneous 
symmetry breaking, similar to the spontaneous magnetization of a ferromagnet which breaks 
the symmetry of the underlying Hamiltonian, implies the existence of massless Goldstone 
bosons which are the pions. The small pion mass of 140 MeV ensures that the concept of 
chiral symmetry persists as a fundamental feature of low energy hadron physics. 
To leading order in density the in-medium scalar condensate drops linearly with the nuclear 
density with the proportionality given by the pion-nucleon sigma term <jn- The model 
independent low density behavior can be derived from the Hellmann-Feynman theorem [l|. 
To determine the corrections to this leading density dependence a large variety of models 
has been exploited. These were e.g. the Nambu-Jona-Lasinio (NJL) model [l|, Q, [J, various 



versions of the linear sigma model [4 
or recently the Polyakov-NJL model 
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5j, the Quark Meson Coupling model (QMC) [6| 



The Hellmann-Feynman theorem relates the in- 
medium scalar quark condensate with the quark mass derivative of the total energy density. 
The latter quantity can also be calculated within hadron effective field theory , such 

as au type models [lOj and, more microscopically within nuclear many-body theory In a 
relativistic framework the latter leads to the relativistic Brueckner (DBHF) approach [ll, 
[l2,[l^] where the nucleon-nucleon interaction is given by boson-exchange models While 
DBHF allows a quite reliable determination of the nuclear equation of state (EOS) up to at 
least two times nuclear density as a conservative estimate - recent DBHF calculations [l5] 
have been shown to be consistent with astrophysical accelerator based constraints concerning 
their high density behavior - the unknown quark mass dependence of the mesonic 
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couplings and masses introduces large errors when conclusions on the in-medium condensate 



are drawn 



Most hadronic models do, however, not respect chiral symmetry. A more systematic and 
direct connection to QCD is provided by chiral effective field theory (EFT). Up to now the 
two-nucleon system has been considered up to next-to-next-to-next-to-leading order (N 3 LO) 

nnrx 

in chiral perturbation theory [17|, [18|, |19j. In such approaches the nucleon-nucleon (NN) 
potential consists of one-, two- and three-pion exchanges and contact interactions which 
account for the short-range contributions. The advantage of such approaches is the system- 
atic expansion of the iViV interaction in terms of chiral power counting. The expansion is 
performed in powers of (Q/A X Y where Q is the generic low momentum scale given by the nu- 
cleon three-momentum, or the four-momenta of virtual pions or a pion mass. A x ~ 47r/ w ~ 1 
GeV is the chiral symmetry breaking scale. Moreover, chiral EFT has a well defined quark 
mass dependence. For the NN interaction the chiral limit of vanis hing current quark and 
pion masses has been evaluated up to next-to-leading order (NLO) [20j, l2l(] and applied to 
the deuteron problem 22]. 

Exploiting the Hellmann-Feynman theorem chiral EFT allows thus a precise determination 
of the scalar condensate at NLO. 

Within finite density QCD sum rules the scalar condensate determines automatically the 
shift of the effective nucleon mass M*. In the sum rule approach scalar and vector fields 
arise naturally from the structure of the quark propagator which is proportional to the 
corresponding condensates. The quark correlation function can be expressed to leading 
order in terms of the scalar condensate (qq) PB and the vector condensate (q'q) PB which is 
introduced by the breaking of Lorentz invariance due to the presence of the medium. The 
identification of the correlation function with the in-medium nucleon propagator of a dressed 
quasi-particle leads to scalar and vector self-energies S s and S D which are of the same order 



in the condensates _ 
The scalar and vector self-energies S s and E Q are on the other hand key quantities of each 
relativistic hadronic field theory. They determine the nuclear EOS, the effective nucleon 
mass M* = M + S s , the single particle potential U cent ~ S D + S s ~ —50 MeV and the spin- 
orbit potential Us.o. oc (E Q — E s ) ~ 750 MeV in finite nuclei. Relativistic phenomenology 



implies large fields of opposite sign E s ~ —350 MeV and E Q ~ +300 MeV 



same is obtained by relativistic many-body theory, i.e. DBHF 11 



12 
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24|]. The 



151 ] . One has, 
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however, to keep in mind that the individual scalar and vector components are interpolating 
fields which do not directly manifest in experimentally accesable observables. Nevertheless, 
when modern nucleon-nucleon interactions are mapped on a Lorentz covariant operator ba- 
sis (using projection techniques as briefly described in Section IHip they reveal large scalar 



and vector fields of comparable size in nuclear matter as a model independent fact [25|, |26 |. 
This holds not only for manifestly covariantly formulated interactions, such as relativis- 
tic One-Boson- Exchange models (Bonn, CD-Bonn, Nijmegen), but also for non-relativistic 
interactions (Argonne Vi 8 , Reid93, Idaho N 3 LO, Vj ow k) as soon as the symmetries of the 
lorentz group are restored. From the analysis of the chiral EFT (Idaho N 3 LO) potential 
I7I . Q we found that these fields are generated mainly by contact terms which occur at 
next-to-leading order in the chiral expansion and which generate the short-range spin-orbit 



interaction 



26]. 



At moderate nuclear densities the N 3 LO scalar and vector fields were found to be in almost 
perfect agreement with the prediction from leading order QCD sum rules |26| ]. The coin- 
cidence of the nucleon mass shifts obtained from QCD sum rules and relativistic nuclear 
phenomenology has been stressed in many works. However, this agreement is not yet fully 
understood. As pointed out in [27] a naive direct dependence of the nucleon mass M* on the 
scalar quark condensate leads to contradictions with chiral power counting. Long-distance 
physics from virtual pions, i.e., the non-analytic term in the expansion of gives a sizable 
contribution to the modification of the in-medium quark condensate. Such contributions 
are, however, found to play only a minor role for the reduction of the nucleon mass. 
In the present work we calculate both quantities, the scalar quark condensate and the 
effective nucleon mass M* from the same chiral effective interaction. The condensate is 
determined making use of the Hellmann-Feynman theorem and the fact that, at least up to 
NLO, the quark mass dependence of the potential is known from its analytic formulation in 
the chiral limit. The effective mass, on the other hand, can be determined in Hartree-Fock 



approximation making use of projection techniques on a relativistic operator basis 25|, |26 |. 
Since the quark mass dependence is known for the chiral effective interaction up to NLO 
this allows furthermore to investigate possible changes of the properties of symmetric nuclear 
matter up to this order. Therefore the large attractive scalar and repulsive vector self-energy 
components as well as the nuclear equation of state are studied in the chiral limit. This has 
also been done using projection techniques on a relativistic operator basis in Hartree-Fock 
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approximation. Naturally one assumes that the magnitude of these fields persists even in 
the chiral limit since hadronic properties are not expected to change dramatically in the 
case of massless quarks or pions. The reason for this is that the expansion of the nuclear 
force in the context of chiral perturbation theory is well defined for small quark masses and 
should still be valid in the limit m q — * which is equivalent to — > 0. 
In the first part of the paper a short description of the chiral effective NN interaction is 
presented. This is followed by the formalism for the relativistic self-energy components in 
Hartree-Fock approximation. The determination of the relativistic self-energy is based on 
a projection technique which allows to transform any two-body potential amplitude onto a 
covariant operator basis or, in other words, to restore the symmetries of the Lorentz group. 
In the following section the implications for the self-energy components and the equation 
of state for nuclear matter are shown when going to the chiral limit. In the last section 
we first present the prediction of the scalar quark condensate in matter derived with the 
help of the Hellmann-Feynman theorem. This is done within two different approximations, 
namely Hartree-Fock and the Brueckner ladder approximation. The latter allows to study 
the influence of short range correlations. Finally the connection between the effective nucleon 
mass and the scalar quark condensate is discussed. 



II. THE EFT NUCLEON-NUCLEON INTERACTION 

Both, the investigation of the structure of the self-energy in nuclear matter and the resulting 
EOS in the chiral limit (m w — > 0) as well as the determination of the scalar condensate in 
matter depend crucially on the exact knowledge of the implicit and explicit current quark 
mass dependence of the nuclear force on which the many-body approaches are based on. We 
are now in the situation to make use of the chiral NN interaction derived in [3] which allows 
an extrapolation in the pion mass where the quark mass dependence is known analytically 
up to NLO . 

The chiral EFT potential consists of one-, two- and three-pion exchanges and regularizing 
contact interactions describing the short-range contributions. The chiral expansion of the 
NN interaction is performed by organizing the contributions in terms of powers of (Q/A X ) I/ 
where Q is the generic low momentum scale given by the nucleon three-momentum, or the 
four-momenta of virtual pions or a pion mass and A x ~ Anf^ ~ 1 GeV is the chiral symmetry 
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breaking scale. 

In H the light quark masS dependence of the nuclear force has Been derived up to next- 
to- leading (NLO) order in the framework of a modified Weinberg power counting, i.e., addi- 
tionally to the one-pion exchange (OPE) potential and contact terms the leading two-pion 
exchange (TPE) has been considered. 

The explicit form of the chiral effective NN potential Vnlo we use is given by 

v NLO = V OPE + V TPE + V cwt , (1) 

where 

yoez = J Si L +2A _ M Jl8 ) r , . T2 (hJMil 9) , (2) 



t/TPE Ti ■ T 2 J r , v 

1 = -384^^ (g) 



4m 2 (5<7 A - 4g 2 A - 1) + ^(23^ - 10^ - 1) + 



4777.2 + <f 2 . 

fgMn^^-lOrj 2 -lU (3) 



m 



- ^L(g) + In jcTi • qa 2 • q - q 2 B\ • cr 2 J 



64tt 2 /, 

ycont = Cs + ^ . ^ + ^2 ^ _ J^gy*^ - 5fl£ + 2) In ^ (4) 

+ ^ (^ - 64^ (16/ '^ - + 2) ln ^) ^ • ^ 

+ dg 2 + C 2 fc 2 + (C 3 g 2 + cS\ti ■ B 2 ) 

+ iC 5 ffl + °" 2 . (£ x g) + C y 6 (g-o ! i)(g-'?2) + C , 7(A;-o ; i)(A;-o ; 2), 

with £4 and the physical values of the nucleon axial coupling and the pion decay constant, 
respectively. Because at NLO any shift in g& and f n for a different value of m n in the TPE 
is a N 4 LO effect for the TPE the physical values gA = 1-26 and f n = 92.4 MeV are used. 
The value of the pion mass is indicated by compared to the physical one denoted by m^. 
L(q) is given by 



lAml + q 2 a/4777 2 + q 2 + \q I 
Hg) = L m = ^L_ lB ^_ . (5) 

A represents the relative shift in the ratio gj^j 'fn compared to its physical value since they 
show an implicit dependence on the pion mass 

(W/x)m, W 



The low-energy constants (LECs) C$,t and £>s,t appear at LO and are related to the Cs,t 
from [19|] via 

C s ,t = C s ,t + mlD s ,T ■ (7) 



The LECs D$,t have not been fixed by experiment till now. In Ref. 20] natural values have 
been assumed for these constants 

D s ,t = , 2 ' 2 , where a s ,r ~ 1 and A x ~ 1 GeV. (8) 

J* x 

The LECs d\6,di8 and U are related to pion-nucleon interactions. We take U = 4.3 which is 
fixed with relatively small error bars. The LECs d\Q,dis are not yet uniquely fixed, i.e. there 
exists a certain range of possible values fixed from different observables. The implications 
on the results induced by the uncertainties of the LECs will be discussed later in greater 
detail. 



III. CO VARIANT REPRESENTATION AND THE NUCLEON SELF-ENERGY 

In the following the mean field in nuclear matter is determined by calculating the relativistic 
self-energy E in Hartree-Fock approximation at tree level. It provides therefore a qualitative 
rather than a quantitative description of the nuclear many-body problem. Nevertheless the 
scale of these fields is already set at tree level. Although essential for nuclear binding and 
saturation, higher order correlations, in particular short-range correlations, change the size 



of the fields by less than 25%, as has been estimated in [26| comparing HF to a full self- 
consistent relativistic DBHF calculation. However, for nuclear binding these deviations are 
essential. To meet the empirical saturation point of nuclear matter one has to introduce a 
self-consistent scheme and to account for short-range correlations. Such calculations have 
to be based on the in-medium T-matrix (or G-matrix) rather than the bare potential V. 
G-matrix correlations will be discussed in connection with the chiral condensate in the next 
section. 

Here we shortly sketch the formalism for the evaluation of the self-energy. More details can 



be found in Refs. [26|] and [131 ]. 
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The self-energy for a nucleon with four-momentum k follows from the interaction V by 
integrating over the occupied states q in the Fermi sea 



T, a p(k, k F ) 



d 4 q 



G° a (q) V A (k,q) aa .,p T 



(9) 



The evaluation of the Hartree integral is sufficient because the matrix elements are fully anti- 
symmetrized containing direct (Hartree) and exchange (Fock) contributions. To determine 
the self-energy only positive-energy states are taken into account as done in the standard 
DBHF approach. The Dirac propagator 



G D (q) = \i + M}2iri5(q 2 - M 2 )Q(q )Q(k F - |q| 



(10) 



describes the on-shell propagation of a nucleon with momentum q and energy E n = 
y/q 2 + M 2 inside the Fermi sea. Due to the G functions in the propagator only positive 
energy nucleons are allowed in the intermediate scattering states which prevents the oc- 
curence of divergent contributions coming from negative energy states. 
Based on symmetry considerations in isospin saturated symmetric nuclear matter the self- 
energy can be written as a sum of a scalar S s , a time- like vector S Q and a spatial vector 
part S v . Thus, in nuclear matter rest frame the Dirac structure of the self-energy has the 
simple form 



E(k, k F ) = E 8 (fc, fcp) - 70 E (fc, k F )+j-k E v (fc, fc F ) 



(11) 



The Lorentz components of the self-energy operator (Q are then given by 13| 

d 3 q M 



(27T)3E q 

k F d 3 q 

9s 



4gs ~9s + a 9a - 



(A^ - q^ 2 



AM 2 



-9 



PV 



(2vr) 3 
k F d 3 q k • q 
(2vr) 3 |k| 2 E q 



29A + T n AM 2 g ™ 



(12) 



- 2 9A + 



k z - q») 2 



q z AM 2 

To evaluate the self-energy operator in the nuclear matter rest frame the two-body inter- 
action matrix V determined in the two-particle centre-of-mass (cm.) frame and usually 
given in the | JLS) basis, has to be represented covariantly by Dirac operators and Lorentz 
invariant amplitudes 28(], a procedure which can be applied to any two-body amplitude. 
This is also the most convenient way to Lorentz-transform the interaction matrix from one 



frame into another 



29|]. 



Naturally a fully relativistic treatment invokes the excitation of anti-nucleons. However, 
standard NN potentials (even OBE type potentials such Bonn, CD-Bonn or Nijmegen) are 
restricted to the positive energy sector and neglect the explicit coupling to anti-nucleons. 
As a consequence one has to work in a subspace of the full Dirac space. This shortcoming 



can be avoided using fully covariant potentials which exp 



30, |31[. The present investigations and those in Refs. 



25 



icitely include anti-nucleon states 



261 ] have, however, been restricted 



;o "standard" potentials based on the no sea approximation. Similarly, the EFT potentials 
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18] applied here and in our previous investigations do not explicitly include anti-nucleons, 
in contrast to covariant approaches which require renormalization procedures to restore chiral 
power counting [32] . 

Working in the positive energy subspace, symmetry arguments and the restriction to on-shell 
scattering allow the two-body matrix elements to be are expanded in terms of five Lorentz 
invariants. A possible choice of a set of five linearly independent operators are the scalar, 
vector, tensor, axial-vector and pseudo-scalar Fermi covariants T m = {S, V, T,P,A} with 

S = 1®1, V = 7 ^®7 M , T = o-^®<v, P = 7 5® 75 , A = 757^ ® 7 57m . (13) 



The choice of the operator basis is not unique. In 



13] it has been shown that the so-called 



complete pv representation is an appropriate choice where the set of covariants originally 



proposed by Tjon and Wallace 



is given by 



T m = {S, -S, (A - A), PV, -PV} . (14) 

PV and PV are the direct and exchange pseudo-vector covariants, analogous to the pseudo- 
scalar covariant P, however, with 75 replaced by — ^)/2M^^. A PV vertex suppresses 
the coupling to negative states and it is consistent with soft pion theorems based on chiral 
symmetry considerations. 

Thus the on-shell (|q| = |q'|) scattering matrix is given by 

V r/ (q',q) = E^(lql^) r - . (15) 

m 

where 9 is the cm. scattering angle and / = 0, 1 the isospin channel. For the Hartree-Fock 
self-energy it is sufficient to consider 6 = when anti-symmetrized matrix elements are used 
since 9 = ir contains then only redundant information. The transformation of the Born 
amplitudes from an angular-momentum basis onto the covariant basis ffl5l) is now standard 



and runs over the following steps: \LSJ) — > partial wave helicity states — > plane wave helicity 



33] . The 



states — > covariant basis. The first two transformation can e.g. be found in Refs. 
last step has to be performed numerically by matrix inversion 13j, |29 |. 
The expressions for the NLO chiral potential given in Eqs. (j2j)-(j3j) are independent of the 
nucleon mass. Nevertheless the nucleon mass appears in the expressions for the calculation of 
the relativistic mean fields, Eq. (|T2|) . as well as in the procedure of projecting the two-body 
amplitudes on the covariant basis. Moreover, the quark mass dependence of the nucleon 
mass is also the leading term which determines the quark mass dependence of the EOS. The 
EOS, i.e. the energy per particle E/A is given by 

d 3 k 



E/A = — / 

PR JF 



p B Jf 2tt 3 

with the single particle potential U s-p _(k, k F ) defined through the fields 

MS, „ E v k 2 



(16) 



U s . p .{k, hp) = — S s - -i^- 



Vk 2 + M 2 S ° 



(17) 



Vk 2 + M 2 

The integration in Eq. ( |T6l) runs over the Fermi sea F and we account for the full momentum 
dependence of the self-energy components S s , S , S v . 

Within the framework of chiral EFT the physical (vacuum) nucleon mass M can be expressed 
as 

M = M + a N (18) 

where Mo is the value of the nucleon mass in the chiral limit. The sigma term <jn represents 
the contribution from explicit chiral symmetry breaking to the nucleon mass and determines 
the quark mass dependence of the nucleon mass 



a N 



E 

q=u,c 



dM 
dm n 



m 



dM 



(19) 



which, through m 2 ~ m q translates into a dependence on the pion mass. The chiral limit of 
the nucleon mass and of the sigma term, respectively, has been evaluated up to NNLO 341 ] . 
where the corrections to the NLO dependence were, however, found to be small. 
In order to be consistent with the NN interaction we account in the following for the pion 
mass dependence of the nucleon mass at NLO (expressions given in 34j) when the self-energy 
components and the EOS are studied in the chiral limit. 
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FIG. 1: (Color online) The tree level scalar (dashed lines) and vector (full lines) self-energy com- 
ponents in matter at kp = 1.35 fm _1 are shown for different values of the pion mass m T . The 
NLO one-pion exchange (left panel) and two-pion exchange (right panel) are shown. The one-pion 
exchange is obtained with d%e = —1.23 and dis = —0.97. 

IV. NUCLEAR MATTER IN THE CHIRAL LIMIT 



The analysis of the EFT potential performed in [26] revealed that the scalar/vector fields 
are generated by contact terms which occur at NLO order in the chiral expansion. These 
are four-nucleon contact terms with two derivatives which generate the short-range spin- 
orbit interaction. The strength of the corresponding low energy constants, in particular 
those connected to the spin-orbit force, is dictated by P-wave NN scattering data. Pion 
dynamics as well as LO and N 3 LO contacts provide only corrections to the fields generated 
by the NLO contact terms. Thus one could expect that the quark mass dependence of the 
fields is mainly determined by the quark mass dependence of the contact terms which is 
moderate. However, before coming to the full self-energy, we consider the contributions 
from one-pion (OPE) and two-pion-exchange (TPE) separately. Fig. [1] shows the scalar 



11 



S s and vector S Q self-energy components from the next-to-leading order OPE and TPE 
contributions at a Fermi momentum of fcp = 1.35 fm _1 which corresponds to a nuclear 
density of ps = 0.166 fm -3 . As a well known result, at the physical pion mass the scalar 
and vector self-energy components from the pseudo-vector OPE are of the same sign and of 
moderate strength. This is also true in the chiral limit. The self-energy components S s , S Q 
approach a constant value of about 30 MeV. For the not yet uniquely fixed LECs entering 
into the expression for the renormalized OPE the values die — —1.23 and dig = —0.97 have 
been taken [20j. The uncertainty due to these LECs di6,i8 does not significantly affect the 



scalar and vector fields. The same is true for the corresponding EOS (see below). This is, 
however, not the case what concerns the scalar quark condensate as discussed in detail in 
the following chapter. 

The scalar S s and vector E Q self-energy components generated by the TPE are already small 
for the physical case (m n = 138 MeV) and are further reduced by m 2, 5 MeV in the chiral 
limit. At zero momentum both components almost vanish and show a slight increase with 
increasing momentum. As for the OPE the fields are repulsive and approach a constant 
value in the chiral limit. 

Next we will consider the role of the contact terms. The contact terms connected to the 
LECs C1...7, Eq. (Jlj), do not depend on the pion mass at NLO. Since the magnitude of the 
scalar and vector self-energy components is mainly set by contact interactions connected 
to the spin-orbit force where the strength is proportional to the LEC C5 in Eq. (Ill) the 
modification of the fields in the chiral limit can in total be expected to be moderate. The 
pion mass dependent part of the contact interactions, i.e., the first to lines in Eq. (jlj) provides 
only small contributions. 

The uncertainties due to the not yet uniquely fixed LECs rfi6,is entering the renormalized 
OPE (J2j) do not strongly affect the self-energy components. However, a second source of 
uncertainty appears in the the part of the contact interactions connected to the not known 
LECs Ds,t which depend on the pion mass. In Ref. 20j this range of uncertainty has 
been explored through an independent variation of the parameters ots,T in Eq. (jSj) in the 



range of —3.0 < o>st < 3.0. In Ref. 



20j this rather wide variation of the LECs D$,t was 



motivated by a wide range of possible parameter sets of NLO LECs fitted with different 
cut-off combinations. However, in the present case - using the Idaho chiral potential - we 
are restricted to one parameter set with a general cut-off of A = 500 MeV. A variation of the 
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k [MeV] k p [fm" 1 ] 

FIG. 2: (Color online) Left: momentum dependence of the tree level scalar and vector self-energy 
components in nuclear matter at kp = 1.35 fm -1 evaluated for different values of the pion mass 
m n . Right: tree level scalar and vector self-energy components in nuclear matter as a function of 
the Fermi momentum k-p for different values of the pion mass m n . 

LECs over a wide range is therefore likely to overestimate the uncertainty originating from 
the LECs D,s,t- Therefore we restrict the present discussion to values as,T ~ 1- Results 
turned out to be stable against a variation of as and ay in the same direction, i.e. small 
deviations from combinations of as,T where both parameters are close to each other do 
practically not change the results. 

In Fig. [2] the full tree-level self-energy components are now shown as a function of the 
momentum k. An approximately vanishing pion mass, i.e. m n = 2 MeV and = 5 MeV, 
leads to a small reduction of the repulsive vector field (~ 30 MeV) and of the attractive 
scalar field (~ 50 MeV), respectively. The same can be seen on the right hand side in Fig. [2] 
where the fields are shown as a function of the Fermi momentum kp. An increase of the 
pion mass to = 200 MeV leads to an opposite behavior. 

In summary, a careful analysis of the chiral EFT NN interaction leads to large scalar and 
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vector fields which essentially maintain their strength in the chiral limit, however, with the 
tendency of a slight decrease of absolute size. 

In this context it may be interesting to compare this behavior with the naive assumption 
of a dropping a meson mass within the framework of Quantum-Hadron-Dynamics (QHD) 
jlo| . In this case the scalar and vector fields are inverse proportional to the masses of the a 
and uj mesons 

£ s = -4rPs > z = +^jpb (20) 

mi mf, 

(7 U 

where g a and g u are the corresponding meson- nucleon coupling constants and ps ~ Pb is 
the scalar nucleon density. The assumption of dropping a and u meson masses according 
to a naive Brown-Rho scaling [35] together with fairly constant couplings would lead to a 
strong increase of scalar and vector fields in size. Chiral EFT predicts the opposite behavior, 
namely slightly decreasing fields. Interpreting this result in terms of the simple QHD picture 
means that the ratio of coupling functions and meson masses in Eq. ( 1201) has to stay fairly 
constant. Assuming dropping meson masses the coupling functions should show the same 
density dependence. Such a scenario is not completely unrealistic since in the framework 
of density dependent relativistic mean field theory [36( where density dependent meson 
coupling functions gi u {pB) are derived from the Brueckner G-matrix [37|, [38] or fitted to 
finite nuclei such a behavior is usually obtained. 

It is clear that a rather small reduction of the scalar and vector fields has only moderate 
consequences for the EOS in the chiral limit. How the change in the fields affects the 
nuclear EOS is depicted in Fig. [31 However, before coming to the pion mass dependence 
we shortly discuss the tree level EOS derived from the chiral EFT potential. The EOSs for 
isospin symmetric nuclear matter at the various orders of the potential are shown on the 
left hand side of Fig. [3j There appear large jumps in the EOS when going from LO up to 
N 3 LO. As discussed in [2jj| the contact terms which generate the large attractive/repulsive 
scalar/vector potentials arise at NLO. Contributions from higher order provide corrections 
to these potentials. These are moderate on the scale of the fields of several hundred MeV 
magnitude. However, due to the subtle cancellation between scalar attraction and vector 
repulsion such corrections may be large on the scale of the binding energy, i.e. several tenth 
of MeV. This behavior is exactly reflected in the EOSs shown in Fig. [3] which jump from 
unbound at LO to over-bound at NLO and N 2 LO to loosely bound at N 3 LO. 
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FIG. 3: (Color online) Hartree-Fock calculation of the nuclear equation of state, i.e. the energy 
per particle E/A as a function of the Fermi momentum kp. On the left hand side the tree level 
results are shown order by order up to N 3 LO. On the right hand side the pion mass dependence of 
the EOS at NLO is shown again, i.e. the result for the physical case of m n = 138 MeV is compared 
to the case of m n = 2 MeV, m I = 5 MeV and m n = 200 MeV. 

The present tree level calculation is of course not a realistic microscopic nuclear matter 
calculation which would on the one hand require to apply the N 3 LO force, and secondly, 
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to perform a self-consistent summation of the Brueckner ladder diagrams 
difference between a tree-level and a full relativistic Brueckner calculation has been discussed 
in 26( . However, the NLO tree level calculation allows a consistent investigation of the chiral 



limit at the order at which the pion mass dependence of the chiral NN potential has been 
derived. 

This is done in the left panel of Fig. [3] which compares the NLO EOS at the physical pion 
mass = 138 MeV to that at m n = 200 MeV, m„ = 2 MeV and m n = 5 MeV, respectively. 
Doing so, one observes a slight softening of the repulsive NLO EOS when the pion mass is 
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increased from its physical value to = 200 MeV. In the chiral limit the scalar attraction is 
slightly stronger reduced than the vector repulsion (see Fig. [2]). The EOS becomes therefore 
more repulsive. One can expect that the observed effect survives also when a full Brueckner 
ladder is summed. Since the pion mass dependence of the contact terms is rather weak, 
Brueckner short range correlations can not be expected to change the results dramatically. 
As already pointed out, the main source of uncertainty arises from the unknown LECs Ds,t 
entering the renormalized contact forces. Concerning the tree level EOS, we did, however, 
not find any tendency of a qualitative change even for large, probably unrealistic variations 
of the dimensionless coefficients ols,t- In a full Brueckner calculation the attraction may be 
increased since iterated OPE and TPE is known to be quenched by Pauli blocking of the 
intermediate states in the Bethe-Goldstone, respectively the Bethe-Salpeter equation. 
Therefore the conclusion is, that the magnitude of the large scalar and vector fields in matter 
persists in the chiral limit and that the physics of infinite nuclear matter is similar to that 
of a vanishing pion mass. This finding is in agreement with the investigations made in 41]. 
In this work nuclear matter was analyzed taking basically the chiral limit of the OPE and 
assuming that the short range and the intermediate range part of any iViV potential is not 
affected. 



V. IN-MEDIUM SCALAR CONDENSATE 

The spontaneous breakdown of chiral symmetry involves a qualitative rearrangement of 
the QCD ground state, due to the appearance of scalar quark- ant iquark pairs. The cor- 
responding non- vanishing ground-state expectation value (qq), denoted as the scalar quark 
condensate, is an order parameter of spontaneous chiral symmetry breaking. Any reduction 
of the scalar density of quarks in matter can therefore be interpreted as a sign of partial 
restoration of chiral symmetry. 

A. Determination from Hellmann-Feynman theorem 



The vacuum value of the lowest- dimensional quark condensate is about [421 ] 



(qq) ~ -(225 ±25 MeV) 3 . (21) 
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The density dependence of the chiral condensate (qq) can be extracted exploiting the 
Hellmann-Feynman theorem with respect to the symmetry breaking current quark mass 
term of the QCD Hamiltonian. We consider isospin symmetric matter making thereby use 
of the isospin symmetry of the condensates ((qq) = (uu) ~ (dd)). Defining qq = ^(uu + dd) 
and m q = \(m u + mA the quark mass term is given by 2m q qq. Isospin-breaking terms are 
neglected. With the help of the Hellmann-Feynman theorem one obtains the in-medium 
quark condensate by determining the energy density £ of nuclear matter 

dS 



2m q ((qq) PB - (qq) ) = m q 



dm n 



(22) 



The derivative is taken at fixed density. The energy density of nuclear matter is given by 

(23) 



S = Mp B + ^1 PB 



A 

where the second term of £ is the energy per particle E/A (times the baryon density), i.e. the 
contributions from the nucleon kinetic energy and nucleon-nucleon interactions. Inserting 
the energy density 8 into Eq. (T22]) and using the Gell-Mann, Oakes, Renner (GOR) relation 



2m q (qq) = -m n f n 

with the definition of the pion-nucleon term from Eq. (TT9l) 

dM a N 



dm q m q 



one obtains 



(qq) 



Pb 



1 



PB 



a N + m, 



d E 



9 dm q A 



(24) 



(25) 



(26) 



(QQ)o ' ™ 2 Jl 

The derivative of the energy per particle with respect to the quark mass can be re-expressed 
using the chain rule 

(QQ)pb _ i 



Pb 



d(E/A) dM d(E/A)dm 7T 
a N + m q \ ' ' — — + m q v 1 1 



dM dm n 



(qq) mlfl 
The derivative of the pion mass using Eq. (|24|) is given by 

dm^ rrijr 



dm n dm q 



(27) 



dm q 2m q 



valid to leading order in chiral perturbation theory. Introducing 



a N 



(28) 



(29) 
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one finally obtains 



x 



d(E/A) d(E/A) m % 



dM dm w 2a n 



(30) 



The first term on the r.h.s. of Eq. fl30l) which reduces the condensate in matter is model 
independent and of first order in the nuclear density [H, 



23]. Inserting the empirical value 
of a N = (45 ± 7) MeV for the sigma term [43j and taking m n = 138 MeV and f n = 92.4 



MeV one finds in Table [U that the in-medium scalar condensate is to leading order in 
density approximately | smaller than its vacuum value at nuclear saturation density. In the 
following a value of po — 0.173 fm~ 3 , corresponding to a Fermi momentum of hp = 1.37 
fm _1 , is chosen as the standard value for the nuclear saturation density. 
From Fig. H] one sees that to leading order a complete restoration of chiral symmetry would 
already occur at ps ~ 2.7po- Such a scenario is unrealistic and contradictory to the knowl- 



edge from heavy ion reactions |44j and astrophysics, e.g. neutron stars 16|. Hence, one 
has to account for higher order corrections in density coming from the d{E / A) / dm q term 
in Eq. (126]) . One might estimate this correction to be small due to the binding energy of 
E/A ~ —16 MeV, which is almost two orders of magnitude smaller than the nucleon mass 
contributing dominantly to the energy density in Eq. ( |23l) . Nevertheless, since the quark 
mass derivative of the interaction energy is the relevant quantity, it is by far not obvious 
that higher order corrections are negligible. 

Thus, a reliable extraction of the density dependent scalar condensate (qq) PB requires both, 
a sophisticated nuclear matter calculation and the exact knowledge of the current quark 
mass dependence of all model parameters entering into the energy density. 
Previous estimates of the scalar condensate based on sophisticated ab-initio many-body 

n q 

approaches [8|, |9J suffered from this problem. The relativistic Brueckner approach chosen 



m 



9| provides a reliable description of nuclear matter bulk properties. Such calculations 



are based on realistic NN potentials, e.g. one-boson exchange potentials [14|. However, the 
current quark mass dependences of the model parameters, i.e. meson masses and coupling 
constants, are unknown to large extent and have therefore either been roughly estimated or 
even been neglected 

In js] it was found that the largest uncertainty in the calculation of the in-medium quark 
condensate (qq) PB arises due to the unknown quark mass dependence of the scalar isoscalar 
a meson exchange which parameterizes effectively correlated two-pion exchange. 
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In the present work the determination of the energy per particle E/A will be determined 
within Hartree-Fock (HF), Eq. ffl6|) . and in a second step within the Brueckner-Hartree- 
Fock (BHF) approximation. Since the NN interaction is thereby based on chiral EFT, 
Eqs. (j2j)-(jlj), where the complete pion mass dependence is known up to NLO [20| we are 
free of uncertainties concerning unknown quark mass derivatives. Remaining ambiguities 
when applying the Hellman-Feynman theorem are only due to the not yet uniquely fixed 
LECs Ds,t i n the NLO contact terms, see Eq. (j3J), and the LECs G?i6,i8 showing up in the 
OPE exchange. The uncertainties coming from these LECs will be discussed. 
In Hartree-Fock approximation nuclear matter is normally unbound, in particular when 
high precision OBE type potentials are applied. The situation turns out to be qualitatively 
different for low momentum interactions like V\ ow k an d Idaho N 3 LO (see Fig. |3D where the 
hard core is strongly suppressed by high momentum cut-offs. In this case isospin saturated 
nuclear matter collapses and Brueckner ladder correlations do not improve on 45|, |46j. In 



this case the repulsion generated by three-body forces which appear at N 2 LO turned out to 
be essential to stabilize nuclear matter and to obtain reasonable saturation properties 46]. 
The inclusion of three-body forces is beyond the scope of the present work, in particular 
what concerns the pion mass dependence, but for a quantitative determination of the EFT 
EOS they should be included. 

However, for a reliable estimate of the in-medium condensate the role of NN correlations, in 
particular short-range and tensor correlations, has to be considerd. One might assume that 
NN correlations influence the result for the condensate, in particular at higher densities. 
In order to estimate their importance the self-consistent iteration scheme of BHF theory is 
applied. 

The central equation of the BHF approximation is the Bethe-Goldstone equation 

G(u) = V + V— %- G(u). (31) 

V is the bare interaction and Q the Pauli operator which prevents from scattering into 
occupied intermediate states below the Fermi momentum kp. The starting energy is denoted 
by uj. The operator Hq defines the energy spectrum of the intermediate two-particle state 
where we use the so called continuous choice. Therefore the single-particle energies for 
particles as well as for hole states above the Fermi surface are calculated from the kinetic 
energy and a mean field part which has to be determined from the Q matrix self-consistently. 
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The single-particle energies are then given by 

e B a HF = e a +J2 (an\G{u = e» HF + e« HF )\a»). (32) 

fi<F 

These single-particle energies are then parameterized in terms of an effective mass and a 
constant potential, i.e. €k ~ k 2 /2M* + U. 

In the present work we apply the non-relativistic approach since it is not possible to use the 



chiral NN potential in a relativistic BHF calculation (DHBF) 13|, ll5j where one accounts 
in addition for the dressing of the potential matrix elements V i— ► V*. The latter requires, 
however, a definite relativistic structure of the interaction, like for covariant OBE-type 
potentials. Nevertheless, differences between a relativistic and a non-relativistic treatment 
should be moderate concerning the derivative of the EOS with respect to the current quark 
mass. 

The prediction of the in- medium scalar condensate in both approaches, i.e., HF and BHF 
are shown in Fig. HI As expected, deviations from the leading order result due to NN 
interactions and nucleon kinetic energy, Eq. fl30l) . increase with density. For both approaches, 
HF and BHF, the additional contributions lead to a weaker reduction of the in-medium quark 
condensate. Especially in the case of the BHF calculation the leading order prediction 
provides a very good description of the quark condensate up to a density of 0.8pB- 
At nuclear saturation density po the reduction of the in-medium quark condensate is about 
3% (BHF) and 12% (HF) smaller compared to leading order. Deviations are, however, 
growing with density, where at ps ~ 2 4- 3po the quark condensate is reduced to ~ 35% 
(HF) and « 30% (BHF). Naturally the BHF approach is more reliable in this density region. 
However, in summary effects from short-range NN correlations and the quenching of OPE 
and TPE due to Pauli blocking, both present in BHF, have only minor implications for the 
condensate as can be seen from the comparison to the HF result. 

The uncertainty due to the not yet uniquely fixed LECs di6,is m the renormalized 
OPE, Eq. ([2]) which was already mentioned in the context of the EOS in the chiral limit, 
enters also into the determination of the scalar condensate. However, now this uncertainty is 
much more severe. The light shaded bands in Fig. H] indicate the range of possible variations: 
The LEC dig is extracted from the Goldberger-Treiman discrepancy. We take the three em- 
pirically found values given in 20| extracted from three different phase shift, d ls = —0.84 
GeV" 



43|, dig = -0.97 GeV" 2 [48j and d w = -1.54 GeV" 2 [49J, respectively. Furthermore 
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FIG. 4: In-medium scalar quark condensate as a function of density with = 45 MeV obtained in 
various approximations: A Hartree-Fock calculation of {qq) PB /{qq)o is shown compared to M*/M 
where M* = M + S s is evaluated at tree level in Hartree-Fock approximation (left). A full 
calculation based on the Brueckner-Hartree-Fock approach is shown on the right. The light shaded 
bands indicate the uncertainty of the corresponding result varying the LECs d±Q from d\Q = —0.91 
GeV" 2 to die = -1-76 GeV" 2 and d 18 from d 18 = -0.84 GeV~ 2 to d 18 = -1.54 GeV" 2 . Dashed 
line: model-independent leading order result. 
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we vary the LEC die in the range from <i 16 = —0.91 to d 16 = —1.76 as done in 
upper bound of the shaded band corresponds to d w = —1.76 and di 8 = —0.84 whereas the 
lower bound corresponds to d\§ = —0.91 and dig = —1.54. These uncertainties are also 
given in Table [H The HF (dash-dotted line) and BHF (solid line) mean values are obtained 
with di 6 — —1-23 and dig = —0.97. 

Comparing with previous approaches performed in a similar spirit 8|, |9J we find generally a 
stronger reduction of the scalar condensate. In [9J calculations were based on the relativistic 
DBHF approach and an one-boson-exchange potential (Bonn A) has been used. In this 
approach an unexpected increase of the in-medium scalar condensate at densities above 
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TABLE I: Predictions of (qq) PB /(qq)o obtained with the Hellmann-Feynman theorem in diverse 
approximations compared with M* /M for three different values of the nucleon density pb- 



Pb/po 


Hellmann-Feynman 

Leading order HF 


theorem 

BHF 


M*/M 
in HF 


0.5 


0.815 


0.828 ± 0.002 


0.812 ±0.001 


0.759 


1.0 


0.630 


0.677 ±0.010 


0.641 ± 0.004 


0.546 


1.5 


0.445 


0.550 ± 0.020 


0.510 ±0.014 


0.354 



p < 2.5pb has been found. The same tendency, i.e. an increasing quark condensate at high 

n 

density has been observed in [8|. In [9] it was assumed that this increase is caused by a 
breakdown of the underlying assumptions related to the current quark mass dependences of 
the model parameters, i.e. meson masses and coupling constants. The authors concluded 
that the use of not chirally invariant NN potentials may lead to wrong predictions in a 
density region where chiral restoration is expected to occur. 

As already mentioned we do not face such problems since the chirally invariant EFT inter- 
action used in the present work has a well denned quark mass dependence. The only source 
of uncertainty arise due to the LECs die ,is which are not yet uniquely fixed and the un- 
known LECs D$,t entering the short-range part, i.e., the contact force which could provide 
substantial corrections to the scalar quark condensate. Nevertheless, the same argument 
given in the previous section kept us from showing a wide undefined variation. Therefore, 
we restrict our calculation again to the case of ols,t ~ 1- However, as for the EOS, the 
prediction of the quark condensate is not considerably altered varying cts,T f° r combinations 
of ols,t where both parameters are close to each other. The contributions which change 
the condensate originate then mainly from TPE and renormalized contact forces. The LECs 
C1...7 in Eq. (jlj) do not depend on the pion mass after renormalization and the related contact 
terms do therefore not contribute to the change of the in-medium quark condensate. 
Nevertheless, considering the possible band of variation due to the LECs di6,i8> both, the HF 
and BHF calculations shown in Fig. H] do not indicate saturation or even an increase of the 
condensate in the considered density range up to 2>ps- Extrapolating the BHF prediction 
to high densities a complete restoration of chiral symmetry, i.e., a vanishing scalar quark 
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condensate is not likely to happen below 4po, even if one takes the range of uncertainty from 
not yet exactly known LECs into account. 

The first determination of the in-medium quark condensate adopting the Hellmann-Feynman 
theorem has been carried out by Cohen et al. [1], based on the tt-N Fock term. There the 
condensate was found to be reduced to a value of 0.694 at nuclear saturation density po 
and 0.58 at 1.5po which is in fair agreement with 0.677 ± 0.01 and 0.550 ± 0.020 obtained 
in the present HF calculation (third column of Table [I]). This agreement, is, however, 
somewhat accidental since we find that TPE and contact interactions (which both have not 
been included in reduce the contribution from OPE by ~ 50%. Moreover, in [lj the 
dg-xN /dm q dependence has been neglected and a different value for g wN has been used. Both 
calculations are, however, comparable in the sense that NN correlations are neglected and 
they are of the same order in the density. Moreover, short-range physics due to contact 
terms, which have been neglected in [l| are found to provide only moderate corrections as 
can be seen from Fig. [5l 

In Fig. [5] the derivatives d(E/A)/dm 7T from Eq. ( 1301 as a function of the density are shown 
with respect to the different contributions coming from pion dynamics, Eqs. ([2]) and ([3]), and 
from the contact terms (V c on, Eq. (j4j), respectively. The dashed-dotted line indicates the 
derivative d{E / A) / DM . The contribution coming from OPE is negative and considerably 
larger compared to those from TPE and contact interactions. In general the contributions 
from pion dynamics, i.e., OPE and TPE are smaller in BHF due to quenching effects. In 
the case of the contact interactions one has to keep in mind that the short-range terms 
rr4 In rfi^ show up due to TPE and the renormalization of the leading order contact terms 
by pion loops. 

The contribution from nucleon interactions is getting substantially more important with 
increasing density compared to the contribution from the nucleon kinetic energy (dot-dashed 
line) which is of order 0(p 5 / 3 ). Thus the nuclear interaction provides important corrections 
to the Fermi gas approximation usually made in QCD sum rule approaches. 
It turns out that the calculations are highly sensitive to the pion mass dependence of the 
pion nucleon coupling constant which has been often neglected in earlier works. This 
fact can also be seen from the relatively large bands of uncertainty in Fig. H] since the 
corresponding LECs enter into the relative shift of gj\/F n , Eq. ([6]), which is connected to 
g-xN via the Goldberger-Treiman relation g^/M = g^/F^. 
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FIG. 5: The derivatives d(E/A)/dm 7T as a function of density with respect to the full NLO cal- 
culation as well as the separate contributions, i.e. OPE, TPE and contact interactions are shown. 
Additionally the dash-dotted line denotes d(E/A)/dM for the full NLO calculation. 

We conclude that the contributions from nucleon interactions to the change of the scalar 
condensate in matter are mainly due to low-momentum virtual pions. In contrast to the 
scalar/vector fields which are generated by NLO contact interactions 25| the contact terms 
and short-range correlations, i.e., the short-distance physics, seem to play a minor role for 
the change of the in-medium quark condensate. Nevertheless, for a fully reliable prediction 
of the in-medium quark condensate the little known LECs D$,t entering the NLO contact 
interactions, Eq. (j4j), have to be fixed with better precision. 



B. Effective nucleon mass 



QCD in-medium sum rules relate the scalar and vector in-medium condensates (qq) PB and 
(q Jf q) PB to the isoscalar scalar and vector self-energies of a nucleon in matter. Thus the model 
independent leading order result which should be valid at low density determines the density 
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dependence of the effective nucleon mass M* = M + S s within the in-medium QCD sum rule 
approach. The scalar and vector fields S s and S Q arise naturally from the structure of the 
quark propagator which is proportional to the corresponding in-medium quark condensate. 
The quark correlation function which follows from the operator product expansion can be 
written to leading order in terms of scalar {qq) PB and vector condensates (q^q) PB - In contrast 
to the scalar condensate to leading order the vector condensate is exactly known. It is given 
by the quark density in the nuclear matter rest-frame since the baryon current is conserved, 



i.e. (q^q) PN = 3/2p B . In [26| we compared the NLO EFT vector self-energy to the leading 
order sum rule result. As in the case of the scalar self-energy deviations were found to be 
small at moderate densities. The next order in the operator product expansion involves 
four-quark operators and combinations of quark and gluon fields which are often neglected 
due to their highly non-trivial structure. Attempts to fix the density dependence of higher 



order contributions in the operator product expansion have e.g. been performed in [501 ] . 
Here we restrict the present discussion to the scalar field S s which follows automatically 
identifying the correlation function with the in-medium nucleon propagator of a dressed 
quasi-particle [l[ Q] 

8vr 2 r . ,_ . . 8tt 2 a N 



A 2 [(ft)?,, - (<i<i)o\ = - tf- — — ps- ( 33 ) 

A B A B m u + m d 

The expression on the right is usually obtained using GOR and the model independent 
term for the scalar condensate in matter which depends linearly on the nucleon density p B 
in Eq. (1261) . Consequently the expression is of leading order in density where ps is the scalar 
nucleon density. The Borel mass scale A B ~ ~ 1 GeV is the generic low energy scale 
of QCD which separates the non-perturbative from the perturbative regime. It coincides 
with the chiral symmetry breaking scale A x of ChPT. Applying Ioffe's formula [5l| for the 



nucleon mass M ~ — j£-(qq) and the GOR relation one finally obtains 52] 



S.W = ~^f 2 PB , (34) 

IT J IT 

where the difference between the scalar and vector density, p$ and p B , can be neglected at 
low densities k 2 F <C M 2 . For the ratio M*/M follows then 

M*_M + £ S a N 

~~M~ ~ ~ </f B - (35) 

7T J TV 

Naturally, Eq. (|35l) represents the model independent leading order prediction for 
(qq) PB /(qq)oi i- e - ; the leading order term in Eq. ff26l) . Equalizing the two quantities, i.e., 
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scalar condensate and effective nucleon mass, is, however, not as straightforward as Eq. (|35|) 



would suggest 
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271 ] . As already stated in Ref. [27j a direct dependence of the nucleon 



mass M* on the in-medium condensate contradicts chiral power counting. Moreover, the in- 
medium quark condensate contains contributions from low-momentum virtual pions, which 
do not contribute to the properties of the nucleon in matter. 

With the present formalism at hand, we are able to perform a consistent comparison of 
the in-medium scalar condensate, derived from the Hellmann-Feynman theorem, and the 
effective nucleon mass where both are obtained from the same chiral EFT interaction and 
at the same order. 

In Fig.H]also the ratio M* /M is shown with the effective nucleon mass given by M* = M+E s . 
The scalar field S s is determined from the chiral EFT potential at NLO in HF approxima- 
tion, making use of projection techniques on a relativistic operator basis as described in 
Section [In] As one can see also from Table (U at saturation density the effective mass M* or 
the ratio M*/M, respectively, is reduced to a value of about 0.546 and is decreasing approx- 
imately linear up to 2po- The reduction of the effective mass M* is about w 13% larger at 
Po than that of the scalar condensate (HF), see also Tabled! At 1.5 po the difference is about 
20 ± 2%. Thus the approximation of Eq. (|34"|) does not hold. By a naive comparison of the 
two quantities the in-medium condensate may contribute at the utmost by about ~ 80% to 
the change of the the nucleon mass in matter at 1.5 po- 

As already mentioned, the higher order contributions in Eq. ( J301) from the nucleon inter- 
action are mainly due to OPE and TPE, i.e., low-momentum virtual pions give the main 
contribution to the change of the scalar quark condensate. The appearance of the large 
scalar field S s which enters the effective nucleon mass M * = M + E s originates on the other 
hand from NLO contact interactions (to be more precise from the part which is connected 
to the spin-orbit force), i.e., it is driven by short distance physics [26]. Low-momentum pion 
dynamics is negligible concerning the appearance of the large scalar and vector fields S s and 
S at the considered order (NLO). The present investigations confirm thus the considerations 
of Ref. [27] which were based on a chiral expansion of the sigma term. 
In summary, a direct dependence of the properties of the nucleon mass on the in-medium 
scalar quark condensate as suggested by Eq. (|35|) can be ruled out. 
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VI. SUMMARY 



We investigated nuclear bulk properties in the chiral limit — > 0. This concerns both, 
the EOS as well as scalar and vector self-energy fields in matter. The large relativistic self- 
energy components are obtained by restoring the symmetries of the Lorentz group of the 
nuclear interaction. The essential ingredient on which the present investigations are based is 
the chiral EFT nuclear force where the implicit and explicit pion mass dependence is known 
analytically up to NLO and allows a well defined extrapolation in the pion mass 20|. 
We found that nuclear bulk properties remain fairly stable in the chiral limit. This is true for 
the EOS as well as the magnitude of the scalar and vector mean fields. Both quantities are 
mainly affected by the pion mass dependence of the OPE and TPE. In such investigations 
a source of uncertainty remains due to the not completely constrained LECs Dg^r which 
appear in the renormalized contact forces. Nevertheless a qualitative change concerning the 
properties of the EOS in the chiral limit is not likely to happen. 

Furthermore we have calculated the density dependence of the chiral order parameter or 
scalar quark condensate in nuclear matter making use of the Hellmann-Feynman theorem 
which relates the scalar quark condensate with the current quark mass derivative of the 
nuclear energy density. As above, the energy density was calculated in Hartree-Fock ap- 
proximation. However, to be more realistic and to include also short range correlations we 
applied also the Brueckner-Hartree-Fock approximation. Since the quark mass dependence 
of the chiral NN interaction is known up to NLO this approach is free from any ambiguities 
which arise concerning the analytic and chiral structure of the potential. The quark mass 
dependence of the pion-nucleon coupling constant g n ^ has thereby been taken into account 
and was found to be important. Uncertainties due to unknown low-energy constants entering 
the pion-nucleon coupling constant g n N do not change the results qualitatively. 
Higher order corrections from the nucleon kinetic and interaction energy become significantly 
more important above saturation density when compared to the model independent leading 
order prediction for the scalar quark condensate. They lead in general to a weaker reduction 
of the in-medium quark condensate and do not indicate a complete restoration of chiral 
symmetry in the density range where hadronic models are reliable. Since Hartree-Fock and 
Brueckner-Hartree-Fock provide quantitatively comparable results one can conclude that 
short-range correlations and quenching effects, both present in Brueckner theory, have only 
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minor implications for the density dependence of the quark condensate. The substantial 
contributions from nucleon interactions are due to low-momentum virtual pions, i.e., OPE 
and TPE. Short-distance physics in terms of contact terms and short-range correlations have 
no important impact on the in-medium properties of the quark condensate. 
This present formalism allows also to perform a consistent comparison of the in-medium 
scalar condensate, derived directly from the Hellmann-Feynman theorem, and the effective 
nucleon mass M* = M + S s where the scalar self-energy S s enters. For the first time both 
quantities were derived from the same chiral EFT interaction and at the same order. In 
general the effective nucleon mass calculated in the many-body approach is smaller (« 10% 
at p ) then the model independent leading order prediction which is used in the QCD sum 
rule approach. 

Moreover, in the present investigations it turns out, that the reduction of the two quantities, 
namely the in-medium condensate and the in-medium nucleon mass, are of different physical 
origin. While the latter is dominantly generated by short distance physics in terms of NLO 
contact interactions 26( virtual low-momentum pions provide the essential contributions 
responsible for the change of the in-medium scalar quark condensate. 
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